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This article reports on our microscopic investigations of the edge of the fractional quantum Hall 
state at filling factor v = 1/3. We show that the interaction dependence of the wave function is 
well described in an approximation that includes mixing with higher composite-fermion Landau 
levels in the lowest order. We then proceed to calculate the equal time edge Green function, which 
provides evidence that the Luttinger exponent characterizing the decay of the Green function at 
long distances is interaction dependent. The relevance of this result to tunneling experiments is 
discussed. 
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There has been interest in the physics of edge states 
in fractional quantum Hall effect (FQHE) jjj for several 
reasons. On the quantum Hall plateaus, the Hall cur- 
rent flows along the edges (at least for sufficiently small 
currents) because that is the only place where gapless 
excitations are present. Presence of extended edge states 
is crucial for the QHE, as is clear in the explanation of 
the physics of the plateaus || . The physics of the edge 
states is also directly relevant to experiments studying 
tunneling into edge states H||. 

Theoretically, the electrons at the edge of a FQHE sys- 
tem constitute an example of a one dimensional liquid, 
the long-distance, low-energy physics of which is generi- 
cally described by the Luttinger model. One of the most 
interesting aspects of the edge Luttinger liquid was the 
assertion, based on general arguments, that the quanti- 
zation of the Hall conductance fixes the parameter char- 
acterizing its asymptotic behavior ||. Consider, for ex- 
ample, the equal time edge Green function, defined as 
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where is the ground state, \& e and \jd are annihilation 
and creation field operators, and x is the distance along 
the edge. Wen argued that the long-distance behavior 
of the Green function is given by \G e dge{x)\ ~ x~ Q with 
a = 3 for v = 1/3; such behavior was verified || in 
an explicit calculation for Laughlin's wave function for 
bosons at v = 1/2. 

However, the exponent measured in the experiment of 
Grayson et al. ^ is a continuous function of the filling 
factor. More recent experiments [Q show evidence for a 
plateau near v — 1/3 but with a ~ 2.7. This raises the 
question that we wish to explore in this work: Does the 
actual form of the interaction affect the properties of the 
edge Luttinger liquid, and if so, how? There has been 
much work investigating this and related issues [(^-[To||, 
and our goal here is to compute the exponent in a micro- 
scopic approach. 



The microscopic understanding of the FQHE is based 
on certain wave functions , which are known to be 

quite accurate for the bulk, homogeneous FQHE states 
from extensive comparisons against exact results in the 
edge-less spherical geometry. However, the accuracy of 
these wave functions remains relatively untested at the 
edges. The validity of the wave functions in the bulk 
does not necessarily carry over to the edges. In the bulk, 
the ground state is separated from excitations by a gap, 
which renders it robust to the actual form of the inter- 
action. On the other hand, at the edges, the system is 
compressible with gapless excitations, which may make it 
more susceptible to various perturbations and the actual 
form of the inter-electron interaction. 

We will work in the disk geometry, which has been 
employed by several authors in the past [pdl[T^ -p^|. We 
neglect in the following mixing with higher electronic 
Landau levels, as appropriate in the large B limit. We 
also do not include in our calculations any external con- 
finement potential; electrons are confined within a disk 
because of the restriction on the total angular momen- 
tum. (These results can be shown to be directly relevant 
to a weak parabolic confinement, which enters into the 
problem only through a renormalization of the magnetic 
length.) 

We will focus on the filling factor at which the edge 
properties have been most investigated, namely the low- 
est Landau level (LL) filling v = 1/3. The bulk ground 
state here is well approximated by Laughlin's wave func- 
tion Hill: 
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•,j — iyj denotes the position of the jth par- 
ticle as a complex number and = y hc/eB is the mag- 
netic length. In fact, is the exact ground state for a 
short range model potential jL8| denoted by V\. However, 
it is not exact for the Coulomb interaction, and signifi- 
cant deviations exist between this wave function and the 
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exact Coulomb ground state in the disk geometry jL7) . 
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FIG. 1. Schematic depiction of the composite- fermion exci- 
ton in which a composite fermion, shown as an electron bound 
to two flux quanta, is excited to a higher CF-Landau level. 
This state is described by the wave function ^ ( ^~ exc ' ton , 
explained in the text, with m* = 3. The y-axis is the com- 
posite- fermion Landau level index (n*), and the x-axis is the 
effective angular momentum m*. 



To gain an understanding of the physical origin of the 
corrections to $J, 3 , we consider the question of mixing 

of ^5/3 w ith low energy excitations, as might be induced 
by the long range part of the Coulomb interaction. We 
construct excitations using the composite fermion theory 
of the FQHE |Q. In this approach, ^Jyg is interpreted 
as one filled composite-fermion LL, rewritten it as 



*?/3 =*i n^-**) 2 

j<k 
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where $i is the wave function of one filled LL, given by 
*i = A[q (r 1 )'q 1 {r2)...ri N - 1 (r N )] 

= Y[( z 3 - z k) ex Ph772 H N 2 ] > 



j<k 
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where A is the antisymmetrization operator and r\ m (r) is 
the single particle wave function in the lowest LL with an- 
gular momentum m. The Jastrow factor fTj<fe (-^J ~~ Zk ) 2 
attaches to each electron in $i to convert it into a com- 
posite fermion. The excitations of ^5/3 are then images 
of excitations in $i. In particular, a particle hole pair 
excitation in $i produces a particle- hole pair of compos- 
ite fermions, namely a composite-fermion (CF) exciton, 
in the full wave function. We restrict our attention be- 
low to the states containing a single CF exciton [|l9| , and 
furthermore, to those excitons for which the CF-LL in- 
dex changes by unity. (Remember: mixing with higher 
electronic LLs is not considered in this work.) The valid- 
ity of this approximation will be checked below. These 
excitons are related to excitations in which one electron 
in <&i is promoted from the angular momentum m* state 
in the lowest electronic LL to the angular momentum m* 
state in the second electronic LL, producing an excited 
state denoted by 



el — exciton 



Almoin)... rj m . _ i (f m * ) Cm* (r 

m*+l )^]m* +1 y m* +2 

(5) 

where Cm* is the wave function of an electron in the sec- 
ond LL in angular momentum m* state. The excitations 
that do not conserve m* are not relevant because they 
have a different total angular momentum than / 3 and 
therefore do not couple to it. The corresponding CF ex- 
citon at v — 1/3 is shown schematically in Fig. ([!]), with 
the wave function given by: 



CF — exciton 
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where Vlll denotes projection of the wave function into 
the lowest LL. The explicit form of the projected wave 
functions can be written in the standard manner [po[ . 
The largest value of m* , corresponding to the angular 
momentum of the outermost occupied state in <3>i, is N — 
1, giving a total of N CF-exciton states. 



TABLE I. The overlap 



1/31*1)3 
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as a function of N, the number of particles. The 
states are explained in the text. Also given are 
E° =< \E-; /3 |#|*; /3 > / < *? /3 K /3 > and 

E mod =< y™°d\ H \y™od > j < ^od^vaod > ^ q{ 

e 2 /eio, where H is the Coulomb interaction. The quoted 
energy is the total energy of the interacting electron sys- 
tem (without taking into account any neutralizing back- 
ground) . 
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We next proceed to diagonalize the Coulomb Hamilto- 
nian in the restricted Hilbert space of N + 1 states: ^5/3 
and ^-ezciton m * = 0) _^ N _^_ Qne difficulty 

is that these states are not mutually orthogonal. We 
obtain an orthogonal basis following the Gram-Schmid 
procedure, compute the matrix elements of the Coulomb 
Hamiltonian, and then carry out diagonalization in the 
restricted Hilbert space, performing all integrals by the 
Monte Carlo method. While the procedure is rather cum- 
bersome, it can be carried out satisfactorily, and we re- 
fer the reader to literature for further details and also 
for a trick used to evaluate off-diagonal matrix elements 
by Monte Carlo To obtain adequate accuracy, we 

need to perform up to 10 million Monte Carlo iterations 
for each point. The modified composite fermion ground 
state obtained in this manner, i.e., taking into account 
the mixing of the "unperturbed" CF ground state ^fy 3 
with single CF excitons, will be denoted by ^i° 3 - We 
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note that like *&y 3 , ^i/| d also involves only single parti- 
cle states with m < 3(N — 1); the states with non-zero 
occupation of m > 3(N — 1) would contain two or more 
CF-excitons, which we neglect in this work. We have 
studied several systems up to N — 30; the overlap be- 
tween \J r ™| d and *j/ 3 , also evaluated by Monte Carlo, is 
shown in Table I. The density profile of ^J/g is shown 
m Figs. (§) and (|) for several values of N. 

0.5 i . . ■ . 1 



we compare in Fig (H) the modified density profile calcu- 
lated above with the exact Coulomb density profile taken 
from Ref . jlTj . The comparison demonstrates that 4"™^ 

is a significant improvement over ^5/3' an< ^ ^ s actually 
quite close to the exact wave function, thus establishing 
that much of the discrepancy between ^5/3 an< ^ the ex ~ 
act wave function can be accounted for by incorporating 
mixing with the lowest energy single CF-excitons. This, 
to our knowledge, is the first variational improvement of 
the FQHE wave function in the lowest LL space. 



N = 10 



10 




FIG. 2. The density profile of y 3 (dashed line) and the 
modified composite fermion state \E r ™| d described in text 
(dots) for N — 10, 11, and 12 particles. The distance from the 
center of the disk, r, is quoted in units of the magnetic length, 
lo = sj hc/eB. The statistical uncertainty in Monte Carlo is 
smaller than the size of the dots. Also given is the density 
profile for the exact Coulomb ground state (solid line), taken 
from Tsiper and Goldman |Tij |. 

To test the validity of the approximations made above, 
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FIG. 3. The density profiles for ^5/ 3 (dashed line) and 
*™3 d (dots) for N = 30. The solid line is a fit according 



to p(r) = 1/3 + 



■ cos ( 2ty 



XQ — OLX I 



with x = r — ro, 



(1 — px-^-yx'^ ) 

r = 11.407, pi = 0.113, x = 4.02, a = 0.029, (3 = 0.415, 
and 7 = 0.086. (All distances are quoted in units of lo.) 

Having demonstrated that ^73 correctly captures the 
interaction dependent physics at the edge, we proceed to 
calculate the edge Green function. For the disk geometry, 
the equal time edge Green function is defined as 



Gedgeifi) 



(7) 



where $ is the ground state describing a disk of electrons 
centered at the origin, z = R is the positions at the edge 
along the x-axis and z = Re = Re l6 is the position of the 
edge at an angle 9 relative to the x-axis. Following Lee 
and Wen |J , we take R — y/6N Iq at v = 1/3. The power 
law exponent a is defined by the equation appropriate for 
the geometry under consideration ||: 

\G edge (9)\~sm- a (0/2) (8) 

The Green function can be expressed as 

J d 2 z 2 .:d 2 z N y*{Re,z 2 , z N )V(R, z 2 , ...,z N ) 



G edge {9) = 



J d 2 z 1 ...d 2 z N ^*(z 1 , ...,z N )V(zi, z N ) 

(9) 
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We write the numerator as <C 'Sg* where <C^ de- 
notes integral over coordinates Z2---ZN- Using 



< >=< aW-a* + i3?*2* > 



2 < 3M>g* >=< |* e +*| 2 > - < |*| 2 > 



< |*e| 2 > 



2 < 3*2* >=< 1*6 



-i*| 2 > - < |*| 2 > - < |* e | 2 > 



we express the real and imaginary parts of the edge Green 
function in terms of integrals over positive definite, real 
integrands, which can be evaluated efficiently by Monte 
Carlo. 

Fig. (|) shows \G edge (9)\ for both *° /3 and *™/° d for 
N = 30. For the former wave function, we recover a = 3, 
as expected. The deviation of the Green function from 
the power-law behavior at short distances (small 9) seen 
in Fig. is quite similar to that found in Ref. || for 
bosons at v = 1/2. There is also a well defined expo- 
nent for *7/3 d , indicating that the edge continues to be 
described as a Luttinger liquid. However, the Luttinger 
exponent now is a w 2.5. 
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FIG. 4. The edge Geeen function as a function of the dis- 
tance along the edge, parametrized by the angle 6. The dis- 
tance from the center is taken to be R — \/~QN lo, lo being the 
magnetic length. The result is shown for *?/ 3 , which is exact 
for the Vi interaction, and the modified CF wave functions 
for two interactions: the Coulomb interaction (r -1 ) and the 
Yukawa interaction (r _1 e _r//i °). In each case, the straight 
line is a fit with the exponent a shown near the line. 



The change in the exponent is surprisingly large con- 
sidering the smallness of the corrections to *?/ 3 - Several 
points are in order here. It is sometimes stated that the 
arguments leading to a — 3 require a sharp edge, i.e., 
when the density is approximately constant (at its bulk 
value) up to the edge and falls sufficiently sharply to 



zero at the edge. Independent of the question of whether 
and how a sharp edge may be realized, we note that both 
*!jy 3 and *5y3 have rather similar density profiles at the 
edge, indicating that sharp-vs-smooth-edge type consid- 
erations are not likely to have any bearing on the dif- 
ference in the exponents. We have also tested that the 
mixing with the CF excitons does not affect the com- 
pressibility of the state; the gap in the angular momen- 
tum L = 3 (TV — 1) sector of the Hilbert space is only very 
weakly affected by it. In other words, the mixing with 
the second CF-LL merely renormalizes the states of the 
lowest CF-LL, while preserving the single channel char- 
acter of the problem. Finally, while it is impossible to 
prove rigorously that the conclusion based on our finite 
system is valid in the thermodynamic limit, the system 
size (N — 30) is sufficiently large to reproduce the ex- 
pected exponent a = 3 for Laughlin's wave function and 
also to give a well defined exponent for the modified wave 
function. 

One might wonder if a universal exponent is obtained 
at least for all short-ranged interactions. To address this 
issue, we have also considered a short-ranged, Yukawa- 
type interaction V(r) = r _1 exp[— r/lo]- The edge Green 
function is computed as before, also shown in Fig. (Q), 
from which the edge exponent a ~ 2.6 is deduced. This 
suggests that the exponent a = 3 is unique to *j , 3 , which 
has the very special property that all of the zeros of the 
wave function are bound to the electrons; any splitting 
of the zeros away from electrons seems to change the 
exponent. 

We believe that our study provides compelling evi- 
dence that the edge exponent is sensitive to the detailed 
form of the wave function. Because the Hall quantiza- 
tion is not affected by small perturbations in the wave 
function, so long as the gap in the bulk is preserved, one 
must conclude that the Luttinger parameter a is not tied 
to the quantized Hall conductance, and is therefore not 
a fundamental property of the FQHE state. Investiga- 
tions into the origin of the non-universality of the edge 
properties are bound to produce interesting physics. 

Our results are generally consistent with the tunneling 
experiments of Chang and collaborators ||[|]. The the- 
oretical exponent obtained above is smaller than 3.0 for 
both the Coulomb and Yukawa interactions, consistent 
with experiment. (A more accurate estimation of the ob- 
served exponent will require a realistic treatment of the 
confining potential at the edge 0).) Also, unless the ex- 
ponent is completely fixed by the Hall conductance, there 
is no reason to expect it to be precisely constant over the 
range of a Hall plateau. 
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